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Abstract 

We argue that the anomalous magnetic moment of the electron (a e ) can be used to probe 
new physics. We show that the present bound on new-physics contributions to a e is 8 x 10~ 13 , 
but the sensitivity can be improved by about an order of magnitude with new measurements of 
a e and more refined determinations of a in atomic-physics experiments. Tests on new-physics 
effects in a e can play a crucial role in the interpretation of the observed discrepancy in the 
anomalous magnetic moment of the muon (a M ). In a large class of models, new contributions to 
magnetic moments scale with the square of lepton masses and thus the anomaly in a M suggests a 
new-physics effect in a e of (0.7 ±0.2) x 10~ 13 . We also present examples of new-physics theories 
in which this scaling is violated and larger effects in a e are expected. In such models the value of 
a e is correlated with specific predictions for processes with violation of lepton number or lepton 
universality, and with the electric dipole moment of the electron. 
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1 Introduction 



The anomalous magnetic moment of the muon = (g — 2)„/2 is one of the most celebrated 
tests of the Standard Model (SM). Indeed, the high precision of its theoretical and experimental 
determinations makes a powerful test on new physics. The situation has become especially 
intriguing with the ~ 3.5a reported discrepancy between the SM prediction and the experimental 
value [l}|4] 

Aa M = ajf p - 4 M = 2.90(90) x 10" 9 . (1) 

On the other hand, the anomalous magnetic moment of the electron a e has never played a 
role in testing ideas beyond the SM. In fact, it is believed that new-physics contaminations of 
a e are too small to be relevant and, with this assumption, the measurement of a e is employed 
to determine the value of the fine-structure constant a. 

The aim of this paper is to emphasize that the situation has now changed, thanks to 
advancements both on the theoretical and experimental sides. Indeed, the theoretical prediction 
of a e has been refined to an unprecedented accuracy and its experimental value is now known 
with smaller errors. At the same time, good determinations of the fine-structure constant 
have been obtained from atomic physics experiments, providing a value of a that is completely 
independent of the measurements of a e . As a result, a e can now be viewed as a very useful 
probe of physics beyond the SM and the situation is going to become even more promising soon, 
as efforts are underway to reduce significantly both theoretical and experimental errors. The 
most exciting aspect of the story is that a e will soon provide us with a crucial consistency check 
of new-physics interpretations of the alleged discrepancy in a^. Moreover, in certain classes of 
models, if Aa M is caused by new physics, then it is possible to correlate the value of a e with 
various other rare processes violating lepton universality or individual lepton number. 

The paper is organized as follows. In sect. [2] we carefully review the present status of the SM 
prediction of a e and confront it with the experimental measurement. In sect. [3] we discuss future 
improvements in the theoretical and experimental results and show their impact for probing 
new physics. We use the prototype examples of supersymmetry (sect. [4]), of a light pseudoscalar 
(sect. [5]), and of vector-like fermions (sect. [6]) to illustrate generic features of theories beyond 
the SM. Our results are summarized in sect. 



2 Status of the electron g — 2 
2.1 The experimental situation 

The classic series of measurements of the electron and positron anomalous magnetic moments 
carried out at the University of Washington yielded in 1987 the value a^ F = 115 965 218 83 (42) x 
10~ 13 [5j[6]. More recently, a new determination of the electron g— 2 has been performed by 
Gabrielse and his collaborators at Harvard University, with the result [7] 

a^ xp = 115 965 218 07.3 (2.8) x 10~ 13 . (2) 
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The uncertainty of this result, (5a EXP = 2.8 x 1CP 13 , i.e. 0.24 parts in a billion (ppb), is 15 
times smaller than that reported back in 1987. The two measurements differ by 1.8 standard 
deviations. 



2.2 The Standard Model prediction 

The SM prediction a| M is usually split into three parts: QED, electroweak (EW) and hadronic. 
Here we provide a summary of the present status of these contributions. 



2.2.1 QED contribution 

The QED term ag ED arises from the subset of SM diagrams containing only leptons and photons. 
This dimensionless quantity can be cast in the general form [8] 

a?*>=A l+ M (^) + AJ^) + AJ^,^-), (3) 



where m e , and m T are the electron, muon and tau lepton masses. The term A\, arising from 
diagrams containing only photons and electrons, is mass and flavor independent. In contrast, 
the terms A2 and A3, generated by graphs containing also muons and taus, are functions of the 
indicated mass ratios. The muon contribution to ag ED , although suppressed by m^/m 2 ~ 2.34 x 
10~ 5 , is about ten times larger than the experimental uncertainty in eq. ([2]) (the tau contribution, 
suppressed by ml/m 2 ., is of order 10 -14 ). The functions Ai (i = 1, 2, 3) can be expanded as power 
series in a/n and computed order-by-order: A4 = A^ (a /it) + A^ (a/7:) 2 + A^ (a/71") 3 + • • • . 

Only one diagram is involved in the evaluation of the one-loop (first-order in a, second-order 
in the electric charge) contribution; it provides the famous result by Schwinger [9] 

Ci = A? = 1/2. (4) 

Seven two-loop diagrams contribute to the fourth-order coefficient A^\ one to A^ (m e /m^) and 
one to {m e /m T ), while (m e /mn, m e /m T ) = 0. The exact mass-independent coefficient 
has been known for more than fifty years, A^ = 197/144 + vr 2 /12 + 3C(3)/4 - (vr 2 /2)ln2 = 



—0.328 478 965 579 193 78... 10 , where is the Riemann zeta function. The numerical evalu- 
ation of the exact expression for the mass-dependent coefficient A^ |ll|l2| with the latest CO- 
DATA recommended mass ratios m e /m^ = 4.836 331 66(12) x 10~ 3 and m e /m T = 2.875 92(26) x 
10" 4 yields A^ ] (m e /m^) = 5.197386 68 (26) x 10~ 7 and A^\m e /m T ) = 1.837 98 (33) x 10~ 9 
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The tiny errors are due to the uncertainties of the mass ratios. The sum of the above values 
provides the two-loop QED coefficient 

C 2 = Af ] + A^ (me/m^ + A ( 2 4) (m e /m T ) = -0.328478 444002 55(33). (5) 

The standard error bC<i = 3.3 x 10~ 13 leads to a totally negligible O(10~ 18 ) uncertainty in a® ED . 

More than one hundred diagrams are involved in the evaluation of the three-loop (sixth- 
order) QED contribution. The exact result for the coefficient A^\ mainly due to Remiddi and 
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his collaborators [14], yields the numerical value = 1.181241456 587... . The analytic 
calculation of the mass-dependent coefficient A^\r) for arbitrary values of the mass ratio r 



was completed in 1993 by Laporta and Remiddi [15]. The exact formula contains hundreds 
of poly logarithmic functions, including harmonic poly logarithms |16| . Numerically evaluat- 
ing these analytic expressions with the latest CODATA mass ratios given above, we obtain 
Af\m e /m^) = -7.373 94162(27) x 1(T 6 and A^ ] (m e /m T ) = -6.5830(11) x 1(T 8 . The same 
values can be obtained with the simple series expansions of [l7j, thus avoiding the complexities 
of these numerical evaluations. The contribution to a^ ED of these three- loop mass-dependent 
coefficients is —0.9 x 10 -13 , i.e. about a third of the present experimental uncertainty <5a EXP , see 
eq. Q. The contribution of the three-loop diagrams with both muon and tau loop insertions in 
the photon propagator can be calculated numerically from the integral expressions of |18|. We 
get ^ 6) (m e /m M ,m e /m r ) = 1.909 82(34) x 10~ 13 , a totally negligible O(10~ 21 ) contribution to 
aJ ED . Adding up all the above values we obtain the three- loop QED coefficient 

C 3 = 1.181234016 816(11). (6) 

The error 6C3 = 1.1 x 10 -11 leads to a totally negligible O(10 -19 ) uncertainty in Og ED . 

Almost one thousand four-loop diagrams contribute to the mass-independent coefficient 

A\ , and only few of them are known analytically. However, in a formidable effort that has 

(8) I 1 

its origins in the 1960s, Kinoshita and his collaborators calculated A\ numerically 19 . Their 

(s) I 1 

latest result is A y °> = -1.9106(20), from Ref. [20]. In the same, very recent, article they 

also computed the tiny mass-dependent four-loop coefficients A*tp \m e /m,j) = 9.222 (66) x 10~ 4 , 

A^\m e /m T ) = 8.24(12)xl0~ 6 , and Ap \m e /m^m e /m T ) = 7.465 (18) x 10" 7 . Adding up these 

values one gets 

C A = -1.9097(20). (7) 

The error 5C4 = 0.0020, caused by the numerical procedure, leads to an uncertainty of 5.8 x 10 -14 
in a^ ED . Independent work on the computation of the four- loop coefficient is in progress [21] . 
Very recently, Kinoshita and his collaborators completed the heroic calculation of the 12672 



five-loop diagrams contributing to the tenth-order coefficient A^ 10 ^ 20|22 . Their result is A\ 



,(io) 



9.16(58). They also computed the tiny mass- dependent term 4 ("le/m^) = -0.003 82(39). 
As the tenth-order contribution of the r lepton loops can be safely neglected, the sum of these 
two numbers provides the five-loop coefficient 

C 5 = 9.16(58). (8) 

The error 5C§ = 0.58 leads to an uncertainty of 3.9 x 10~ in cig ED . 
The complete QED contribution is given by 



a? ED 



(a) = y jTCi(aM i . (9) 



i=l 



As (a/71") 6 = 1.6 x 10 16 , terms of order i > 6 are assumed to be negligible at present. 
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2.2.2 Electroweak and hadronic contributions 

The determination of the EW contribution has been recently refined in [2j: 

a™ = 0.3854 (42) x 1(T 13 . (10) 



This precise value includes the two-loop contributions first computed in [23 J . 
The hadronic term, 

< AD = 16.82(16) x 10~ 13 , (11) 

is six times larger than the present experimental uncertainty 5a EXP , see eq. (|2]). It is the sum of 
the following three contributions: the leading-order one, 18.66 (11) x 10~ 13 , very recently updated 
in 



24 



(see also [2]), the higher-order vacuum-polarization part, —2.234(14) x 10 13 24 (see 



also (§[25]), and the hadronic light-by-light term, 0.39 (13) x 10~ 13 (2] (see also [26]). 
2.2.3 Standard Model prediction of a e and value of a 

The sum of the QED contribution a^ ED (a) plus the hadronic and weak terms discussed above 
yields the SM prediction of the electron g— 2: 

af(a) = a« ED (a) + a EW + < AD (12) 

(the dependence on a of any contribution other than aJ? ED is negligible). To compare it with 
experiment, we need the value of the fine-structure constant a. The latest determination of a 



by CODATA 13 



a (CODATA) = 1/137.035 999 074 (44) [0.32 ppb] , (13) 

cannot be employed for our purpose, as it is mainly driven by the value obtained equating the 
theoretical SM prediction of the electron g— 2 with its measured value, 

af» = a EXP (14) 



(thus assuming the absence of any significant new- physics contribution) Q Indeed, solving eq. ( 14 ) 
with the experimental value of eq. ([2]) we obtain 

a (g-2) = 1/137.035 999174(34) [0.25 ppb], (15) 

in agreement with Refs. [20] . This is the most precise value of a available today. The difference 



between this number and q(CODATA) in eq. (13) is mainly due to the very recent QED five-loop 



result of Ref. 20 , which is included in our derivation, but not in the CODATA one. 

Clearly, in order to compute a| M (a) and compare it with a EXP we must use a determination 
of a independent of the electron g— 2 . At present, the two most accurate ones are 

a( 133 Cs) = 1/137.0360000(11) [7.7ppb], (16) 

a ( 87 Rb) = 1/137.035 999 049(90) [0.66 ppb]. (17) 

The identical value presently reported by the Particle Data Group (PDG) [27], adopted from CODATA, 
cannot be used either. 
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They differ by less than one standard deviation. The first value was obtained from the ratio 
h/Mc s (h is Planck's constant and Mq s is the mass of the 133 Cs atom), which was determined 
by measuring the atomic recoil frequency shift of photons absorbed or emitted by 133 Cs atoms 



using atom interferometry 28 . The second was deduced from the measurement of the ratio 
/i/Mpb (-^Rb is the mass of the 87 Rb atom) with an experimental scheme that combines atom 



interferometry with Bloch oscillation |29-3lj. The values of a in eqs. ( 16|17 ) were inferred from 



the ratios /i/-^Cs,Rb combining them with the very precisely known Rydberg constant and the 
mass ratios Mc s ^/m e 13 . Given the higher precision of a( 87 Rb) vs. a( 133 Cs) (by more than 



one order of magnitude), the former is the value of a we employ to compute a^ M (a). We note 



that a( 87 Rb) agrees with a(g— 2) in eq. (15) (the difference is 1.3 standard deviations), and its 
uncertainty 5a( s Rb) is larger than Sa(g— 2) just by a factor of 2.7. 

The SM prediction Og M (a), computed with the fine-structure constant value a( 87 Rb) of 



eq. (|17|), is 

af* = 115 965218 17.9 (0.6)(0.4)(0.2)(7.6) x 1(T 13 . (18) 

The first (second) error is determined by the uncertainty of the four(five)-loop QED coefficient, 
the third one is 5a^ AD , and the last (7.60 x 10 -13 ) is caused by the error 5a( 87 Rb). The 
uncertainties of the EW and two/three- loop QED contributions are totally negligible. When 
combined in quadrature, all these uncertainties yield 5a^ M = 7.64 x 10~ 13 . Note that the present 
precision of the SM prediction, which is about three times worse than the experimental one, is 
limited by the uncertainty of the fine-structure constant a( 87 Rb). 

2.2.4 Standard Model vs. measurement 



The SM value in eq. (18) is in good agreement with the experimental one in eq. (|]). They differ 
by 

Aa e = a* xp - af 1 = -10.6 (8.1) x HT 13 , (19) 

i.e. 1.3 standard deviations, thus providing a beautiful test of QED at four-loop level! (The four- 
loop contribution to a^ ED is —5.56 x 10~ n .) Once again, the uncertainty 5Aa e = 8.1 x 10~ 13 is 
dominated by that of the SM prediction, through the error caused by <5a( 87 Rb). 

3 New Physics tests with a e 

New physics effects in the electron g— 2 are usually expected to be comparable with the EW 



contribution, a® w = 0.3854(42) x 10~ 13 , see eq. (10), and therefore much smaller than the 
uncertainty 5Aa e = 8.1 x 10~ 13 reported above. Indeed, as we mentioned earlier, the com- 
monly used CODATA (and PDG) value of a is mainly derived from Og XP under this assumption. 
However, as we will discuss in the next section, the uncertainty in Aa e is expected to be re- 
duced. Then, the anomalous magnetic moment of the electron will provide us with important 
information on new physics effects. 
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3.1 Future improvements in the determination of Aa e 



As we showed in sections 2.2.3 and 2.2.4, the uncertainty 5Aa e = 8.1 x 10 13 is the result of the 
combination, in quadrature, of the following errors, in units of 10~ 13 : 

(0.6) QE D4, (0.4) QED5 , (0.2)had, (7.6)* a , (2.8) 5(I exp. (20) 



(0.7)th 

The first one, 0.6 x 10~ 13 , is caused by the numerical procedure used to determine the four-loop 
QED coefficient. It is already rather small and can be further reduced to 0.1 x 10~ 13 with a large 
scale numerical recalculation |32|. Also the second error, induced by the five- loop QED term 
recently computed by Kinoshita and his collaborators, is already small, and may soon drop to 
less than 0.1 x 10 -13 [32]. As the tiny - but hard to reduce - hadronic uncertainty is 0.16 x 10 -13 
and those of the EW and two/three-loop QED contributions are totally negligible, the overall 
purely theoretical error of Aa e (i.e. the value of 5Aa e obtained setting to zero 5a^ xp and the 
the error induced by 5a) is, at present, 0.7 x 10~ 13 . This is likely to decrease even further, by 
a factor of two or three, in a relatively near future. 

Thanks to these recent theoretical improvements, the precision on Aa e is now limited only 
by the experimental uncertainties 5a^ xp and 5a. At present they affect 5Aa e by 2.8 x 10 -13 
(<5ag XP ) and 7.6 x 10~ 13 (5a). It seems reasonable to expect a reduction of the former error 
to a part in 10~ 13 (or better) in ongoing efforts to improve the measurement of the electron 



(and positron) anomalous magnetic moment 33 , 34 . Work is also in progress for a significant 



reduction of the latter error 31,35 



In conclusion, a determination of Aa e at the level of 10 -13 (or below) is a goal that can 
be achieved not too far in the future with ongoing experimental work. As we will discuss in 
the next section, this will bring a e to play a pivotal role in probing new physics in the leptonic 
sector. 

3.2 General structure of new-physics contributions 

The one-loop SM electroweak contribution to the anomalous magnetic moment of the lepton £ 

(I = e, fx, t) is 

(«f W )u oop = 7^2 (l - I - 2 *w + I sm 4 %)^2x 10- ^f. (21) 

This can be viewed as a benchmark for contributions from new physics at the electroweak scale. 
For the muon, this is about the same size as the observed discrepancy, see eq. Q. 

New physics effects for the leptonic g— 2 can be accounted for by means of the effective 
Lagrangian 

C = e^j- {i R a^A u ,e' L + I L a^A* w l R ) F» v £, l' = e, Mj r , (22) 



which describes dipole transitions (£ — > l ,n f) in the leptonic sector. Starting from eq. (22), we 
can evaluate Aa^ as 

Aa e = 2mj Re(A u ) . (23) 
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Let us consider now some new particles with typical mass A NP and couplings g\ and gf to 
left- and right-handed leptons £, respectively. The one-loop new-physics contribution to the 
amplitude Ayy is then of the form 



At 



1 



(4^A NP )2 
and therefore Aai reads now 

2m? 

Aa e 



(4ttA np ) 2 



{9lk 9vh + 9ik 9e k ) fi(x k ) + — (gf k gPt) f 2 {x k ) 



(Iffil 2 + I4l 2 ) h{x k ) + —Re {gf k gf k *) f 2 (x k ) 

rri£ 



(24) 



(25) 



With f\ t 2 we indicate loop functions which depend on ratios (x k ) of unknown masses of the 
new particles contributing to the amplitude £ — > i'j, and A; is a lepton flavor index. In the 
term proportional to /i, the chiral flip required by the dipole transition occurs through a mass 
insertion in the external lepton line. In the term proportional to f 2 , the mass insertion is in 
the internal line of some new particle, thus explaining the parametric factor v/me- Although 
f 2 must be proportional to the lepton Yukawa coupling, as a consequence of chiral symmetry, 
in practice this term can become very sizeable whenever a new large coupling leads to a chiral 
enhancement. 



In a broad class of theories beyond the SM, g e 
same for any £) and f 2 vanishes, such that 



L,R 



A «4 



m ii 
mi, 



and /i are flavor universal (i.e. are the 



(26) 



We will refer to this case as "naive scaling" (NS). NS applies, for instance, if the new particles 
have an underlying SU(3) flavor symmetry in their mass spectrum and in their couplings to 
leptons (which is the case for gauge interactions). 

An interesting consequence of NS is that an explanation of the muon g— 2 anomaly makes 
definite predictions for new effects in the anomalous magnetic moments of electron and r, 



Aa P 



Aa, 



3 x 10" 9 

Aa M 
3 x 10- 9 



\ 0.7 x 1(T 13 , 
0.8 x 10~ 6 . 



(27) 
(28) 



It is very exciting that the sensitivity in Aa e is not far from what is required to test whether the 
discrepancy in the muon g— 2 also manifests itself in the electron g— 2 under the NS hypothesis. 
Thus, determining Aa e with a precision below 10~ 13 is an important goal with rich physics 
consequences. The measurement of a T may play a similar role if a precision below 10~ 6 will be 



attained (see sect. 3.3). 



Although the NS case is especially simple and common to a large class of new-physics 
interactions, it is by no means the only possibility. Many theories make predictions beyond 
NS, either because of the chirally-enhanced term proportional to f 2 , or because of lepton non- 



universality in the couplings (g e ' ) or the mass spectrum (/i). For instance, in a multi-Higgs 
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doublet model, the couplings g e ' are related to Yukawa couplings and therefore one would 
find the scaling AaiJ Aai j = m 4 ./m 4 .. A different example is offered by supersymmetry with 
non-degenerate sleptons, in which Aai loses its correlation with lepton masses. 

The case of non-naive scaling is interesting from the theoretical point of view because it 
allows for exploration of the structure of lepton symmetries. It is also interesting experimentally 
because it can lead to effects in Aa e testable already with present sensitivity. Moreover, as we 
will show with various examples in the next sections, any lepton non-universality in Aai can be 
related to other experimental observables, offering the possibility of cross-checking new physics 
effects. 

The underlying £ —> £'"f transition described by the effective Lagrangian of eq. (22) can 
generate, in addition to the anomalous magnetic moments Aai, a ^ so lepton flavor violating 
(LFV) processes (£ / £'), such as /i — > ej, and CP violating effects, such as the leptonic electric 
dipole moments (EDMs, de), 



BR(^ 
BR(£^, 



48vr 3 a 



\Ap 



+ \A V e\ 2 ) , 



dt 

e 



mi Im(^4 



it) 



(29) 



Using the general parametrization for Aw of eq. (24) we find 
BR(^ -> £' 7) 



BR(^ 



di 
e 



(4vrA N p) 
v 



9ik 9e*k + 9tk del) h{x k ) + —gfk 9i?kh{xk) 



+ £^ 



Im {<&<]%) h{xk). 



(30) 



(4vrA NP )2 

On general grounds, one would expect that, in concrete NP scenarios, Aai, dg and BR(£ — > i'j), 
are correlated. In practice, their correlations depend on the unknown flavor and CP structure 
of the couplings g L and g R , and thus we cannot draw any firm conclusion. In the following, 
we will point out the general conditions that have to be fulfilled by any NP theory in order to 
account for large effects in Aai while satisfying the constraints from di and BR(£ — > £ >r y). 
Regarding the leptonic EDMs, we find the following model-independent relations 

Aa P 



d e 



d- T 



7 x 10- 14 

Aa M 
3 x I0- 9 

Aa-r 



10" 24 tan< 



e cm , 



2 x 1CT 22 tan< 



4 x 10 



-21 



tan < 



e cm . 



e cm . 



(31) 



where we have defined (j) 
current anomaly and Aa e 



x 10- 7 

= arg(A^). We have normalized Aa^ to the central value of the 
to their values in naive scaling, Aae/ Aa^ = mf/m 2 ^. 
From eq. (31 ) we learn that an explanation of the Aa^ anomaly implies, for a natural CPV 
phase cp^ ~ 0(1); a model-independent upper bound on du < 3 x 10 _22 e cm which is still 
far from the current bound d^ < 10~ 18 e cm, but well within the expected future sensitivity 
dp < 10~ 24 e cm [36]. Therefore, any experimental effort to improve the resolution on d^ would 
be valuable. 
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On the other hand, the electron EDM imposes a bound on the corresponding CPV phase 
(j) e at the level of 10~ 3 , if NS is at work. Such a condition could be realized for instance if (fi e is 
generated radiatively while (j>n arises already at the tree level. Going beyond NS, one could also 
envisage scenarios where the electronic dipoles are suppressed compared to the muonic dipoles 
because of hierarchical couplings ge ,R <C gjl ,R , as it happens for instance in a multi-Higgs doublet 



model where g* ' are related to Yukawa couplings. In general, as shown by eq. (30), the EDMs 
(but not Aau) vanish if g L = g R as it could arise in a left-right symmetric theory. 

Now we discuss the potential constraints on Aai arising from LFV processes. An inspec- 



tion of eq. (30) shows that BR(£ — > £ ,r y) and Aai are generally correlated as BR(£ — > £ /r y) 
(Aa M ) 2 \9ui\ 2 where 9u/ stands for the relevant flavor mixing angle, which is approximately 



9lk9 lk 



'<i' 9 11 9 S k i see ec L- (|30h . In particular we have 



-12 / Aa-n \ i ", ,, 



3xl0- 9 / V 3xl0_5 

2 / n \ 2 



BR(/i ery) « 2 x 10 
BR,^, 7 ) - 4xl0~^^j . (32) 

LFV contributions to Aai can be particularly important as they are typically chirally enhanced 
by m T /me, as we will discuss in the case of supersymmetry. Such a chiral enhancement is not 
effective in r LFV processes and therefore, in this case, it might be possible to keep BR(t — > ij) 
under control, while generating large effects especially for Aa e . 

Before turning our attention to some prototype theories beyond the SM that predict non- 
naive scaling in Aai, we will now briefly summarize the status of the anomalous magnetic 
moment of the r. 

3.3 Status of the r anomalous magnetic moment 

The present experimental sensitivity of the r lepton g— 2 is only 10 -2 . In fact, while the SM 



prediction for a T is precisely known 37 



a s T M = 117721 (5) x 10~ B , (33) 

the very short lifetime of this lepton (2.9 x 10 -13 s) makes it very difficult to determine its 
anomalous magnetic moment by measuring its spin precession in a magnetic field, like in the 
electron and muon g— 2 experiments. Instead, experiments focused on high-precision measure- 
ments of r pair production in various high-energy processes and comparison of the measured 
cross sections with the SM predictions. 

The present PDG limit on the r g— 2 was derived by the DELPHI collaboration from e + e~ — > 
e + e~T + T~ total cross section measurements at LEP2: —0.052 < a® xp < 0.013 at 95% confidence 
level [38] - This reference also quotes the result in the form: 

a^ xp = -0.018(17). (34) 
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Comparing eqs. (33) and ( p4| ) (their difference is roughly one standard deviation), it is clear 
that the sensitivity of the best existing measurements is still more than an order of magnitude 
worse than needed. In 39 , the reanalysis of various measurements of the cross section of the 
process e + e _ —> t + t~ , the transverse r polarization and asymmetry at LEP and SLD, as well 
as of the decay width T(W —> tu t ) at LEP and Tevatron allowed the authors to set a stronger 
model-independent limit on new physics contributions: —0.007 < a^ p < 0.005. 

The possibility to improve these limits is certainly not excluded. Future high-luminosity B 
factories, such as Belle-II [40j or SuperB |!41~|, offer new opportunities to improve the determi- 



nation of the r magnetic properties. The authors of 42 proposed to measure the r anomalous 
magnetic moment form factor at the T resonance with sensitivities down to 10~ 5 or 10~ 6 , and 
work is in progress for the determination of a T via the measurement of radiative leptonic r 
decays [43]. 



4 Super symmetry and a e 

The supersymmetric contribution to ag comes from loops with exchange of chargino/sneutrino 
or neutralino/charged lepton. Therefore, violations of "naive scaling" can arise through sources 
of non-universalities in the slepton mass matrices. Such non-universalities can be realized in two 
ways. 

1. Lepton flavor conserving (LFC) case. The charged slepton mass matrix violates the global 
non-abelian flavor symmetry, but preserves U(l) 3 . This case is characterized by non- 
degenerate sleptons (mg ^ rnp, ^ m?) but vanishing mixing angles because of an exact 
alignment, which ensures that Yukawa couplings and the slepton mass matrix can be 
simultaneously diagonalized in the same basis. 

2. Lepton flavor violating (LFV) case. The slepton mass matrix fully breaks flavor symmetry 
up to U(l) lepton number, generating mixing angles that allow for flavor transitions. 
Lepton flavor violating processes, such as fi — >■ ej, provide stringent constraints on this 
case. However, because of flavor transitions, a e and a M can receive new large contributions 
proportional to m T (from a chiral flip in the internal line of the loop diagram), giving a 
new source of non-naive scaling. 

We will discuss these two cases separately. 



4.1 Lepton flavor conserving case 

In the LFC case, we assume non-degenerate slepton masses for different families (mg ^ 7^ rrif) 
but flavor alignment between lepton and slepton mass matrices. This is reminiscent of the 



well-known alignment mechanism 44 , which was proposed to solve the supersymmetric flavor 
problem in the quark sector by aligning the down quark/squark mass matrices and which might 
arise naturally in the context of abelian flavor models |45|. 
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In this case the supersymmetric contribution to Aa\ ¥0 is given by the following approximate 
expression 



Aaf 



arrin 



Re(^M 2 tan/3) 



+ 



4vr sin 2 9 W mj(M 2 -fi 2 ) 
am 2 Re(^Mitan/3) 



fn(x 2 ,x ll ) - fc(x 2 ,x,j) 



8?r cos 2 w mj(M 2 



am i 



Re(A e Mi - ^Mitan/3) 



47T cos 6w 



9n{xi) 



e,/i,T. 



(35) 



The loop functions f n , f c and g n (with ^1(2) = ^i(2)/ m % anc ^ x m = t J - 2 / m 'j) are given in the 



appendix, and we use standard notation for the supersymmetric parameters. In our numerical 
analysis we use the exact expressions for Ao£ FC in the mass eigenstate basis [46 1. 



In the illustrative case of a single mass scale (Mi 2 
arg(^Mi) = 0, the result simplifies to [^] 



nip) and if arg^Mi^) 



Aaf 



5ct!2 mj 
487r raj 

3 x 10~ 9 



tan /3 



ay ra| 
247r ra| 



At 



m. 



100 GeVV 



tan /3 



+ 0.1 



A ± 
3m; 



(36) 



This shows that the discrepancy in the muon g— 2 can be explained by supersymmetric particles 
with masses around 100 (300) GeV, for tan/3 = 3 (20). On the other hand, the contribution 
proportional to the term Af, which could induce in principle a violation of the NS if A e /A^ 7^ 
m e /m^, is too small to explain the muon g—2 anomaly after imposing the vacuum stability 
bound \A e \/m^<3. 

Assuming that sleptons are the heaviest particles running in the loop, we find 

Aa„ 



Aa P ~ Aa 



m 2 th3, 

A 1 2 2 

m~ 



2 



Aa T ~ Aa 



m 2 rri] 
mf. m 

r 1 



rrv 
1 

m 



3 x 10" 9 
3 x 10- 9 



10 



-13 



10" 



(37) 



For values of Aa^ explaining the muon g—2, non-degenerate sleptons at the level ra^ ~ 3 rag lead 
to Aa e ~ 10 -12 , which is at the limit of present experimental sensitivity. This naive expectation 
is confirmed by a complete numerical analysis, as illustrated in the left panel of fig. [TJ which shows 
Aa e as a function of the normalized selectron/smuon mass splitting X efJi = ra|)/(ra|+ra|). 
The plot has been obtained by scanning the supersymmetric parameters that account for the 
(g — 2) M anomaly at the level of 1 < Aa M x 10 9 < 5 (black points), or 2 < Aa M x 10 9 < 4 (red 
points in the inner region). We have imposed all current bounds on sparticle masses arising 
from direct searches and required M x < 100 GeV, M 2 = 2M 1 , u, < 200 GeV, 3 < tan/3 < 30, 
< 500 GeV. As can be clearly seen, values of Aa e almost up to 10 -12 are reachable for 



m 



e,fi 

highly non-degenerate selectron and smuon masses such that \X e ^\ ~ 1. 

2 The overall sign of the supersymmetric contribution to ai is a free parameter, since it is determined by 
arg(/iMi,2), which is invariant under field phase redefinition. 
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In the right panel of fig. [I] we show Aa T as a function of X^ T = (m~ — m~) / (m~ + m|). The 
plot has been obtained following the same procedure described for the left panel. We observe 
that Aa T can reach the level of 10 -5 if m? <C mn, therefore violating NS, whose prediction is 
shown by the horizontal line in fig. [T] The condition nif <C nrijx is justified in the so-called split 



family models [47,48 , where the first two generations of squarks and leptons are substantially 
heavier than the third generation. 

Before concluding this section, it is worth to comment about the potential constraints arising 



from the leptonic EDMs, on the same line done in the model-independent analysis of sect. 3.2 



From eq. (31) we learn that Aa e 10~ 13 ( 12 ) is allowed only if arg^Mi^) < 10 -3 ^. Moreover, 
since we are dealing now with flavor-blind CPV phases, that is <j> = (f> e = <j)n = (j> T = ^sifiMi^), 
the electron EDM limit d e < 1.5 x 10~ 27 e cm implies the bounds < 3 x 10~ 25 e cm and 
d T < 5 x 10 _24 e cm, if mg ~ run ~ vtif. Yet, values of cL as high as the model-independent 



upper bound d M < 10 e cm of eq. (31) could be reached if mg 3> m^, in which case the 



supersymmetric contributions to Aa e and d e could be made negligible. 

4.1.1 Correlation between a e and violation of lepton universality in LFC 

In supersymmetric theories, NS violations for the leptonic (g — 2)i can arise through sources of 
non-universalities in the slepton masses. In turn, these non-universalities will induce violations 
of lepton flavor universality in low- and high-energy processes such asP-> £v, r —> Pu (where 
P = tt,K), ti — > Ijvv, Z —7- It and W — > Iv through loop effects. Lepton universality has been 
probed at the few per-mill level so far, see table [TJ It is interesting to study the correlation 
between such violations of lepton universality and departures from NS for Aa^. Taking for 
example the process P — > iv, we can define the quantity 

( ^ )eXP = 1 + Ar^ . (38) 

\Rp )sM 

Here (R e ^) SM = T(P -> eu) SM /T(P -> fj,v) SM and (R e ^) EXP = F{P -> ei/) EXP /r(P iw) 

EXP 

so that Ar e J M signals the presence of new physics violating lepton universality. 

Within supersymmetry, in the absence of LFV sources, Arp ^ is induced at the loop level 
by box, wave function renormalization and vertex contributions from sparticle exchange. The 



complete calculation of [i decay in supersymmetry 49 , 50 can be easily applied to meson decays 




and we have included the full result in our numerical analysis. Neglecting box contributions, 
which are suppressed by heavy squark masses, the parametrical structure of Arp^ is 

(39) 

where the term t; 2 /min(m? ^) stems from SU(2) breaking effects. 

Within supersymmetry, such SU(2) breaking sources arise from left-right soft breaking 
terms, from mixing terms in the chargino/neutralino mass matrices, or from D-terms. For 
highly non degenerate sleptons and if min(m| -) ~ v 2 , Arp^ can reach the few per-mill level. 
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Channel 

T(tt -> eu e )/T{-K -)■ fxPfj) 
F(K -> eP e )/r(K -> 
r(iir-^7reP e )/r(iir-^7r^i/ M ) 
r(Z -> e + e")/r(Z -»• 

r(r -> i/ T e i> e )/r(r v T ^ i? M ) 



-0.0045 ± 0.0032 
0.004 ± 0.004 
-0.002 ± 0.004 
-0.0010 ±0.0026 
0.017 ±0.019 



51 



52 



53 



27 



27 



54 



54 



-0.0036 ± 0.0028 



56 



56 



57 



Channel 



T(tt ->■ /i f Ai )/r(r -> 7r i/ T ) 

r(z -> / u+^-)/r(z -»• r+r - ) 

r(/i -> e u e )/r(r v T e v e ) 



0.016 ±0.008 
0.037 ±0.016 
-0.0011 ±0.0034 
-0.060 ±0.021 27 



57 



27 



54 



54 



56 



56 



fr 

-0.0014 ± 0.0044 57 



Channel 



Ar e/T 



T(Z -> e+e-)/r(Z -»• r+r - ) 
r(W -»• eP e )/r(W -»• rP T ) 
->• f M e v e )/T{r ^ v T n Pfj) 



-0.0020 ± 0.0030 
-0.044 ±0.021 



27 



27 



54 



54 



56 



56 



-0.0032 ± 0.0042 [571 



Table 1: Experimental limits on Ar^ 1 ". 
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Given that the NP sensitivity of ir — > £v and K — > tv to the above effects is the same, 
and that their present experimental resolutions are comparable, both tt ^ £v and K — > Iv 
represent useful probes of this scenario. Future experiments at TRIUMF and PSI aim at testing 
lepton universality in it — > £v at the level of < 1 x 10~ 3 and 5 x 1CP 4 , respectively. Both the 
NA62 experiment at CERN and the KLOE-2 experiment will continue improving their sensitivity 
aiming at a test of lepton universality in K — > £v at the few per-mill level. 

In fig.[2j on the left, which has been obtained through the same scanning procedure described 
for fig. jlj we show the quantity Ar^/^, defined in eq. (38), as a function of Aa e . We learn that 



large NS violations typically imply also breaking effects of lepton universality at the few per-mill 
level, which are expected to be within the future experimental reach. 

Similarly, in fig. [2j on the right, we show breaking effects of lepton universality in the [ijr 
sector, accounted for by the quantity Ar^^ which can be constructed for instance combining 
processes like r — > Pv and P — > fj,v (where P = tt,K). Large NS violations, bringing Aa T to 
the level of 10 -5 , are typically correlated with violations of lepton universality in the /i/r sector 
at the per-mill level. These effects will be tested experimentally at the SuperB, which aims at 



LFU tests in r decays well below the per-mill level 58 



4.2 Lepton flavor violating case 

We now consider LFV contributions to a^. The importance of this case lies in the fact that LFV 
can generate contributions to a e and that are chirally enhanced, being proportional to m T . 
Using the mass-insertion approximation, which is valid for near degenerate sleptons with mass 
rrif, we find 



Re (mM^VED 



Aaf v ~— tan/3 K - " M LL/ U Xl ) £ = e,fi. (40) 

We defined 5 l ^ B = (m^)^ /m| (with A, B = L, R), where m\ L and m 2 RR stand for the left- and 
right-handed slepton mass square matrices. The function l n is given in the appendix. In the 
illustrative case where = M\ , we find 



As evident from eq. ( 40 ) , the single power of mi and the flavor mixing angles of the soft sec- 
tor break NS and provide a potentially large chiral enhancement m T /m etfl , which is especially 
important for Aa e . 

4.2.1 Correlation between a e and r — > e~y in LFV 

The case of LFV has to be confronted with the strong constraints from processes like li —> £j^, 
see table [2] The decay rate of li — > £j*y is given by 

BR(£i -> Li) 48vr 3 a 



BR(£ 4 -> tjViVj) G\ 



i^r + i^r • (42) 
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Figure 1: Left: Aa e function of X efl = (■ 
of = (m~ — m~)/(m~ + m~). Black point; 
red points correspond to 2 < Aa^ x 10 9 < 4. 



''1 ~~ m fi)/( m i + Right: Aa T as a function 
satisfy the condition 1 < Aa^ x 10 9 < 5, while 




Figure 2: Left: Ar p as a function of Aa e , where Ar p measures violations of lepton universal- 
ity in r(P — > ev)/Y(P — > jjiv) with P = K,ir. Right: Ar^ 1 " as a function of Aa T where Ar^ T 
measures violations of lepton universality in T(P —> /j,u)/T(t —> Pv) with P = K,n. Black 
points satisfy the condition 1 < Aa^ x 10 9 < 5, while red points correspond to 2 < Aa^ x 10 9 < 4. 
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BR(/i~ - 


->• e _ 7) < 2.4 x 10" 


12 


|5£ 




BR(r- - 


-> < 4.4 x 10" 


-8 


60 


I 61 




BR(r- - 


-»■ e _ 7) < 3.3 x 10" 


8 


60 







Table 2: 90% C.L. limits on radiative LFV decays. 



The supersymmetric contributions to [i &y and r — >■ £"f (£ = e, fi) amplitudes arise at one 
loop level through the exchange of charginos/sneutrinos and neutralinos/sleptons. Although our 
numerical results are based on the exact expressions for A % [ R in the mass eigenstate basis 



62 



in the following we provide their expressions in the mass-insertion approximation 63 

a 8 % l L fiM 2 tan j3 



+ 



47rsin 2 % mj (M$ - pp) 

- ^^Mitan/3 

47T cos z &w m~ 



+ 



a 



RL 



4-7T cos 2 9w m 2 - 



Afi 

mi. 



' Xfi ) 

h n (xi) g n {xi 

J Xfl j 

mj ~ (M 2 - JP) 

2/n(xi) , 



(43) 



a 



R 4-k cos 2 9w 



%M lt an/3(^ + 2 ^i^ 



mj (M 2 -fi 2 ) 



+ 



mi 



Mi 

m ti 



(44) 



In the region of parameter space that maximizes Aa e , corresponding to moderate/large tan/3 
values and a large fi term, BR(r — > £j) can be approximated as 



BR(r -> ii) « 3 x 10" 



tan^\ 2 /2TeV\ 4 ( ^/m i 
30 J V m ( - 



0.5 



+ 



0.5 



(45) 



Combining eq. (45) with eq. (41), one can set an upper bound on Aae as a function of BR(r 



£7). We have studied this correlation numerically and we show in fig. [3] our results for Ao e vs. 
BR(t —> ej). The plot has been obtained by means of a scan with Mi < 2 TeV, M 2 = 2M\, 
Ii < 5 TeV, 3 < tan/3 < 30, mg < 2 TeV, \S R e R \ < 1, \8 T L e L \ < 1. As illustrated in fig. g| the 
experimental limit BR(t — > &y) < 3.3 x 10~ 8 (90% CL) curiously corresponds to a limit on Aa e 
close to 10 -12 , which is just at the edge of present experimental sensitivity. Therefore, in the case 
of LFV, any future positive indication for Aa e will necessarily imply that r — > ej is just beyond 
the present bound. Moreover, even if r — >• e7 will not be detected at the upcoming SuperB and 
Belle II facilities, which are expected to reach the sensitivity of BR(t — > e*y) <3x 10~ 9 , we can 
still expect values for Aa e up to the level of Aa e <2x 10 -13 , which is roughly a factor of three 
above the expectations of the NS scenario. 

However, in contrast to the LFC scenario, LFV is not able to account for the Aa^ anomaly, 



with departure from NS. This is because by combining eq. (|45j) with eq. (41) we find that the 
experimental limit BR(r — > fj/y) < 4.4 x 10~ 8 (90% CL) sets a very strong upper bound on 
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Figure 3: Left: BR(r — > e~f) vs. |Ao e | in the LFV scenario. The horizontal lines show the 
current bound and the expected future experimental sensitivity on BR(r — > ej). The vertical 
line corresponds to the prediction for Aa e assuming a NS setting Aa^ to its central value Aa^ = 
3xl0 -9 . Right: Aa^ vs. |Aa e | in the LFV scenario imposing the bound BR(r — > ej) < 3.3xl0 -8 . 
The black line shows the correlation between Aa^ and Aa e in the case of naive scaling. The 
horizontal dashed (solid) lines show the la (2a) Aa^ anomaly Aa^ ~ (3± 1) x 1CP 9 , see eq. (II h. 




Figure 4: Ar e J^ vs. |Aa e | in the LFV scenario. The vertical line corresponds to the prediction 
for Aa e assuming NS, setting Aa^ equal to its central value Aa M = 3 x 10~ 9 . 
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Aajf v . Moreover, large effects in Aa^ FV require slepton masses roughly at the TeV, in order 
to keep BR(r — > e^f) under control. In such a regime, we cannot employ the contribution to 
(<? — 2)^ from LFC because Aa^ FC is not large enough in order to account for the anomaly, see 
eq. (36). This can be clearly seen in fig. [3] on the right where we show Aa^ vs. |Aa e | in the LFV 



scenario, imposing the bound BR(r -4 &y) < 3.3 x 1CP 8 . The black line shows the correlation 
between Aa^ and Aa e in the case of NS where LFV contributions are absent. In the interesting 
region where 10 -13 < Aa e < 1CV 12 , it turns out that Aa^ is below the level of 1CV 9 . 

Finally, there is a formidable constraint from /j, — > ej, a process that receives contributions 
from the combination of flavor mixing angles 5^ e ~ 5^ T 5 Te . Therefore, in order to fulfill the 
BR( / u — > ej) bound while generating sizable effects for Aa e oc 5 eT 5 Te , we need a strong hierarchy 
such that 5 I Y B <C $ab' ^ n conclusion, a large contribution to Aa e in the LFV case is possible if 
the only sizable mixing is between the first and third generation of sleptons. 

Also in the LFV case a potential constraint can arise from the leptonic EDMs. In partic- 
ular, in order to obtain Aa e ~ 10 _13 ( 12 ) we need a suppression for the relevant CPV phase 
arg(/zMi(5^ R 5J^) < 10 _3 ( 4 ). However, in contrast with the LFC case, we have now flavor- 
dependent CPV phases and therefore the electron EDM does not constrain directly cL-r- Still, 
whenever c^ iT are induced by flavor- dependent phases (coming from LFV sources), powerful 
bounds are obtained by the LFV processes l{ — > lj~f. In particular, we find that the current bound 
BR(t — > fvy) < 4.4 x 10~ 8 sets the upper bounds d^ < 3 x 10~ 23 e cm and d T < 1.5 x 10 _24 e cm. 

4.2.2 Correlation between a e and violation of lepton universality in LFV 

Violations of lepton universality can arise also in the LFV scenario. The quantity that is deter- 
mined experimentally and accounts for deviations from /i-e universality is 

/exp Li r ( P 

The sums extend over all neutrino (or antineutrino) flavors since they cannot be distinguished 
experimentally. This is important for our purposes because in the LFV case one expects new 
contributions to P — > ev T , a process that is detected as breaking of lepton universality, rather 
than violation of lepton flavor. 

One would naively expect that the LFV channels P — > (i ^ k) are suppressed compared 
with the LFC ones (z = k). An interesting exception is provided by charged Higgs mediated 
LFV contributions, which can be sizable for large tan /3 and which are chir ally-enhanced by the 



factor m T /m£ 64 . Indeed, the dominant contribution to Arp^ in the case of mixing between 



the first and third generation sleptons is 



1 + Arf 1 ~ 



2 

m P m r A ll .3 



M 2 H± m ( 



Ap L tan^ f3 



2+ (^)(3 ||A ^ ta " "- (47) 



The coefficients Ap L and App measure the effective couplings H + v e i J eR and H + v T ien respec- 
tively, which are induced at one- loop level by exchange of Bino or Bino-Higgsino and sleptons. 



19 



Since these effective Yukawa interactions are of dimension four, the quantities and A^ are 
not sensitive to the overall soft scale, hence avoiding supersymmetric decoupling. In the region 
of parameter space where Aa e receives large effects, i.e. for fi S> rug ~ Ml, it turns out tliatj^] 



An inspection of eqs. (40 48) reveals that LFV effects contributing to Aa» FV and Art/ M are 



correlated, as shown in fig. [4} The plot has been obtained by means of a scan with 10 < tan j3 < 
50, M H + < 1 TeV, Mi < 2 TeV, M 2 = 2Mi, fi < 5 TeV, m i < 2 TeV, and \5 T A e B \ < 1. It is 
interesting that Aa^ FV and Ar|^ can simultaneously reach experimentally testable values even 
for supersymmetric masses beyond the LHC reach. 

Turning to pion physics, one finds that Arv^ ~ (m^/m|) x Ar^/ M . Therefore Ar^Jsusy is 
negligible after the constraints from Art/ 1 ?, are imposed. 

4.3 Disoriented A-terms 

So far, we have investigated various supersymmetric scenarios and their capability to account 
for the muon g—2 anomaly. In particular, we have analyzed the possibility of breaking the NS 
among different leptonic g—2. As discussed in the previous sections, both model-independently 
as well as in supersymmetric frameworks, the major challenge we have to deal with, when 
generating large effects for the leptonic g—2, is to keep under control other dipole transitions 
like the electron EDM and [i —> ej. Yet, the correlation among Aa^, di and BR(^ — > i'j) depends 
on the unknown flavor and CP structure of new physics. Therefore, from a phenomenological 
perspective, dg and BR(^ — > i'^f) do not provide a direct bound on new contributions to Aa£. 

However, it would be desirable to have a concrete scenario where the above conditions 
are naturally fulfilled. In the following, we point out that this happens in the case of the so- 
called "disoriented A-terms" , invoked in ref. |66| to account for the recently observed direct CP 
violation in charm decays D — > KK, tttt. 

The assumption of disoriented A-terms is that flavor violation is restricted to the trilinear 
terms 



1 f = u,d,£, (49) 



m f 



where 0- are generic mixing angles. This pattern can be obtained when the trilinear terms 
have the same hierarchical pattern as the corresponding Yukawa matrices but they do not 
respect exact proportionality. A natural realization of this ansatz arises in scenarios with partial 



compositeness, as recently pointed out in ref. |67|. Interestingly, the structure of eq. (49) allows 



us to naturally satisfy the very stringent flavor bounds of the down-sector thanks to the smallness 



In particular, starting from the exact expressions for £± RR and A RL in the mass eigenstate basis 65 , which we 
use in our numerical analysis, one can find that A^j ~ ^z^rr an d A^ ~ — j^LL. -^-SYl^r f° r M ^ m e ~ Mi. 
Notice that lm(8'fJ L 5 R e R ) is strongly constrained by the electron EDM. However, sizable contributions to Ar^ 
can still be induced by Ke(SfJ L S R e R ) . 
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of down-type quark masses. Similarly, also the bounds from the lepton sector can be satisfied 
under the (natural) assumption that the unknown leptonic flavor mixing angles are of the form 



m,- 



67 



In the disoriented A-term scenario, the dominant amplitude for fi — > is 

A L = o fn(*l) , 50 

where x\ = Mf/m'~. Assuming that the only possible sources of CP violation arise from A 
terms, the electron EDM d e is generated by the one-loop exchange of Bino and charged sleptons. 
One can find the following approximate expression 

4 = aIm(WS t ) 
e z7r cos^ 0w m 

On the other hand, the g— 2 does not require any source of CP violation and therefore it always 
receives effects also from SU(2) interactions. In particular, the leading effects is 

a mj tan B ... . , . 

Aa e ~ . < : 2 f{x 2 ) , (52) 
7r sin Uw fn 

where we have considered the illustrative case where M 2 = fJ>, so that x 2 = M^/fh? = fj, 2 /m 2 . 
The loop function is such that /'(l) = 1/8 and /'(0) = 1/2. 

Therefore, the leading SUSY contribution for g— 2 is parametrically enhanced relatively to 
the amplitude generating d e and \x — > ej by a factor of ~ tan/3/ tan 2 9w ~ 100 x (tan/3/30). This 
naturally raises the question of whether it is possible to account for the (g — 2)„ anomaly, while 
satisfying in a natural way the constraints from d e and \x — > &y. Indeed, the latter constraints 
generally require slepton masses around the TeV scale, if we assume A ~ rh with 0(1) phases 
and Ofj ~ \fmi/mj. Even with such large slepton masses, it is still possible to account for 
the (g — 2)^ anomaly if gauginos are significantly lighter than sleptons, which is also welcome 
from naturalness arguments. This can be seen observing that if x 2 <C 1, that is if M 2 ,fi -C rh, 
the relevant loop function for Aa^ is enhanced by a factor of four compared to the case where 
M 2 = ii = rh. At the same time, both BR(/x —> cy) and d e tend to decrease when M\ <C rh, 
sec cqs. ( |50|51| ). Clearly, one could equivalently consider a lighter supersymmetric spectrum to 



explain the muon g— 2 and A% terms somewhat smaller than rh. 

In particular, setting rh = \A e \ = 1 TeV, sinc/) J 4 e =l, M 2 = \jl = 2M\ = 0.2 TeV, and 
tan/3 = 30, we find that 



BR(/j -> 67) « 6 x 10~ 13 



TeV ^rrie/m^ 



V m l J 



o s , /'^ 0fi\ /TeV\ 2 
4 x 10~ 28 Im 11 e cm . 



TeV / \ m; / 



Aa M « 1 x 10 



_,, / TeV\ 2 /" tan ,5 



/ V 30 
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Figure 5: Predictions for \i — > ej, Ao M and d e in the disoriented A-term scenario [66] assuming 
9fj = \frrii/ mj . Left: \i — > vs. Aa^. Right: <i e vs. Aa M . 

These estimates are fully confirmed by the numerical analysis shown in fig. [5] which has been 
obtained by means of the following scan: 0.5 < \A e \/m < 2 with sm<f>A e =l, fh < 2 TeV, 
(M 2 ,/x,Mi) < 1 TeV and 10 < tan/3 < 50. 

It is interesting that disoriented A-terms can account for (g—2)^, satisfy the bounds on 
/i — > ej and d e , while giving predictions within experimental reach. However, we expect that 
the electron g— 2 follows NS. A potential source of NS breaking comes from the trilinear terms 
(provided A e /A^ / m e /m^). In practice, as already discussed in previous sections, their effects 
are very small after the vacuum stability bound \A(\jmi < 3 are imposed and therefore the NS 
relations are preserved. 



5 Light (pseudo)scalars and a e 

In this section we will investigate scenarios where the muon g— 2 anomaly is accounted for 
by contributions arising from light (pseudo)scalar particles. Irrespectively of the underlying 
theoretical motivations, such framework is interesting because it typically predicts large and 
very special NS violations. 

We parametrize the Yukawa interactions between the light scalar field (ft and pseudoscalar 
A with leptons I with the following effective Lagrangian 

where and C^ A are arbitrary constants. Although we will not discuss specific models, the field 
A could arise as a pseudo-Goldstone boson of an extended Higgs sector and the field 4> could 
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be a light gauge singlet coupled through a dimension-five effective interaction to the ordinary 
Yukawa terms. 

Very light scalar or pseudoscalar particles with Yukawa-like couplings are generally subject 
to stringent constraints both from low-energy data, such as meson decays, as well as reactor 
and beam dump experiments (for a review see ref. |68|), Most of these bounds disappear for 
Ma ^ 10 GeV, which is the mass regime of interest for us in order to account for the (g—2)^ 
anomaly, as we will show below. For simplicity, we can also consider the case in which <j) and 
A are coupled to leptons, but not to quarks. Because of the smallness of the electron Yukawa 
coupling, this case is even less constrained. 



5.1 One-loop effects 

The one- loop contribution due to (pseudo) scalar particles to (g—2)i is |69 




32tt 2 M^ 



(55) 



where the loop functions Is a are 



A = f dz / { \ He > J i = t dz ~, T~~2~T" > ( 56 ) 

Jo I — z + z z r^ J 1 — z + z^r 1 



A 



with ri A = mf/M^A and with asymptotic limits 



V - I 2r r > 1 j £ _ ( 2r r > 1 ^ 

" ' lnr-| r«l ' A \ -lnr-^ r<l. 



Note that the one-loop pseudoscalar (scalar) effect is unambiguously negative (positive). 
Therefore, at this level, a light pseudoscalar particle cannot explain the muon g—2 anomaly as 
it contributes to Aa^ with the wrong sign. Moreover, for me M^^a we nn d NS, Aa^ oc mj, 
see eqs. ( 55|57 ). On the other hand, in the opposite limit me <C M^a-i ^ a e roughly scales 



with the fourth power of lepton masses (apart from a mild logarithmic dependence) Aae oc 



mgln(MfaA/nT>£): see eq. (57). 



5.2 Two-loop effects 

Since in the me <C Mj, a regime the one loop contributions to Aae are highly suppressed by the 
fourth power of lepton masses, two loop effects might be relevant whenever we can avoid large 
powers of light lepton masses. This is indeed the case for two-loop Barr-Zee type diagrams with 
an effective vertex generated by the exchange of heavy fermions. Here we consider only 
the effect of the r lepton, but top and bottom should also be included if A (or cf) are coupled 
to quarks. Therefore, the two-loop contribution is 



^ )2Ioop " "SvrW^il;^ ( Re (°^*) Ml ~ Re ( C ^ r ) Mj 
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where the loop functions are 



L d> = dz ~r, \ F ln ? — » L A = dz — - - In , 59 

9 Jo z(l-s)-rj r T , J 2(1 



with asymptotic limits 



61nr+13 r vs. i ( lnr+2 vs. i 

^=<! 3 , Z£ = 2 r 3 (60) 

1 ln 2 r + 21nr + 4+ \ r«l \ln 2 r + ^- r«l. 



As shown by eq. (58), two loop effects for Aa£ exhibit NS and can be positive or negative 
depending on the sign of Re(C^C T *). Moreover, the enhancement factor mj/m^ of two loop 
effects relative to one-loop effects, which is particularly important for the case of the electron 
g—2, can easily compensate the additional loop suppression, as we will see later. On the other 
hand, in the case of the t, one-loop effects are typically dominant compared to two-loop effects. 

Hereafter, we will focus on the case of pseudoscalar particles, which is especially interesting. 
Specializing to the mass regime mi <C Ma where the Aa^ anomaly can find an explanation, we 
have the following situation: 1) Aa e is always dominated by two-loop effects, 2) Aa^ receives 
comparable one- and two-loop contributions, and 3) Aa T is always dominated by one-loop effects. 
As a result, we expect significant NS violations that we are going now to study numerically. 

In the left plot of fig. [6j we show the anatomy of the contributions to Aa^ setting Ca = 50: 
the red line corresponds to the magnitude of the negative ("— ") one-loop effects, the green 
line refers to the positive ("+") two-loop effects and the black line stands for the magnitude 
of the total contribution. The most prominent feature emerging from this plot is the different 
decoupling properties of the one- and two- loop effects as one can check directly from eq. ( 60 ) . In 
particular, two-loop effects have a much milder decoupling with Ma compared to one loop-effects. 
This implies, in turn, the existence of a value for Ma where the two effects have comparable 
(and opposite) size. The exact value of Ma where this happens depends on the masses and 
couplings of the particles circulating in the second loop of the Barr-Zee diagram. In our case, an 
almost exact cancellation occurs for Ma ~ 6 GeV and therefore the total contribution to Aa^ 
is positive for Ma > 6 GeV and negative for Ma < 6 GeV. 

In fig. [7] we show the values attained by Aai induced by pseudoscalar effects monitoring, in 
particular, whether NS is at work or not. 

In the left plot of fig. [7j we show Aa^ vs. Aa e for different pseudoscalar masses Ma and 
varying Ca- As we can see, NS (black line) is systematically violated by a large amount and, 
in particular, Aa e always lies above its naive expectation. The actual amount of NS violations 
depends on the degree of cancellation between one- and two-loop effects entering Aa M . Overall, 
in the regions where the Ao^ anomaly is accommodated, Aa e typically exceeds the 10 -13 level, 
thereby providing a splendid opportunity to test the g— 2 anomaly via the electron one. 

In fig. [7j on the right, we show Aa^ vs. Aa r for different pseudoscalar masses Ma and 
varying Ca- Similarly to the case of the electron g— 2, NS (black line) is largely violated and 
Aa T can reach values up to the level of 10 -3 , while explaining the Aa^ anomaly. Such values 
should be well within the experimental resolutions expected at a SuperB factory. 
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Figure 6: The different contributions to Aa M induced by a pseudoscalar particle (with mass Ma) 
as a function of Ma- The red line corresponds to the negative ("— ") one-loop contribution, the 
green line to the positive ("+") two-loop contribution and the black line to the total effect. 




Figure 7: Left: Aa^ vs. Aa e for different values of Ma varying Ca- Right: Aa^ vs. Aa T for 
different values of Ma varying Ca- The naive scaling predictions are given by black lines. 
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Let us discuss now the predictions for the leptonic EDMs induced by light pseudoscalars. At 



one-loop level, the effective Lagrangian of eq. ( 54 ) leads to a real dipole amplitude and therefore 
the EDMs are vanishing at this order. Two loop effects are generally complex and provide very 
important contributions, as already discussed for the case of the leptonic g— 2. Since two-loop 
effects follow NS, it is expected that the electron EDM is the most sensitive probe of this scenario 
among the leptonic EDMs. In particular, from the model-independent expectations of eq. (31), 
we deduce that, in order to accommodate the (g — 2)„ anomaly while satisfying the electron EDM 
bounds, we need the condition Im(C e C^*) < 1CP 3 , where / stands here for the heavy fermion 
running in the second loop. This condition could be naturally satisfied either if Im(C ) = Im(C^) 
or if a dynamical mechanism suppresses both Im(C^) and Im(C^), for instance by a loop factor. 



6 Vector-like fermions and a e 

In the following, we will consider the impact of heavy vector-like fermions on the leptonic (g — 2) 
(see also [70|)- The introduction of heavy vector-like fermions, mixing with the SM fermions, 
can be motivated by the explanation of flavor hierarchies in the SM, as we will discuss shortly. 

Hereafter, we will focus on the leptonic sector only and therefore we introduce SU(2) vector- 
like doublets (Ll © L R ) and singlets (El © E R ) governed by the Lagrangian 

- C = M E E L E R + M L L R L L + m E E L e R + m L L R l L 

+ X LE L L E R H + \ LE L R E L rf + h.c. (61) 

The fields Ll (Lr) and E R (El) have the same (opposite) quantum numbers as the SM fields £l 
and e R , respectively. Note that the SM Yukawas are assumed to be vanishing and they will be 
generated dynamically after electroweak symmetry breaking through the mixing between light 
and heavy fermions, once heavy fermions are integrated out. 

The above Lagrangian is reminiscent of the fermionic sector of composite Higgs models 
which are dual of warped 5d models. In such a scheme, the chiral fermions correspond to 
weakly-coupled elementary fields while the vector-like fermions are composite fields belonging 
to a strongly interacting sector. The Higgs boson also belongs to the strong sector and does not 
couple to elementary fields. 

The mass eigenstates, before electroweak symmetry breaking, are obtained diagonalizing 
the mass mixing in C through the following 2x2 unitary matrices 




cos 9 L sinflA ( e R \ / cosfljj sinfljA / e R 

-sm9 L cosO L ) \L L J ' \E R J \-sm6 R cos6 R J \E R/ 



(62) 



where 



tan#r, = —r- 1 tan#p = — t~ ■ (63) 
Ml M e y ' 

Hereafter we use the notation sl( R ) = sm9 L ( R ) and cl( R \ = cos Oli R )- After performing the 
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rotations of eq. (62), the Lagrangian of eq. (61) becomes 

-C = M' E E L E R + M' L L R L L + \ LE L R E L H^ 

+ A L £ (s l sr Il&r + c L s R L L e R + s L c R £ L E R + c L c R L L E R ) H + h.c. , 



where 



M'r 



M? + 



M' 



Ml + 



L L — V L L ' 1V1 E — y ±V± E 

and therefore we can define the following mass mixing matrix 



m 



E - 



(64) 
(65) 



h ( ^LESLS R V 



M± 



>^lecls r v 




Lr 

M' L 

>^LEV 




(66) 



M' E 



As a result, the fermionic spectrum consists of two heavy fermions with masses approximatively 
given by M' E and M' L and light SM fermions with masses given by 



mm 



^>LE S L S R V + 



(67) 



where i = 1,2,3 is a lepton flavor index. The leptonic g— 2 are generated at the loop-level by 
the exchange of the Higgs, W, Z bosons and heavy fermions. The dominant contribution arises 
from the diagram with an underlying mixing of SU(2) doublets and singlets (since it is chirally 
enhanced by a factor of vjmi). 



Using the mass relation of eq. (67), neglecting terms of order 0(v ' /Mf E ) and performing 



an explicit loop calculation one can find that the dominant contribution is given by |70 



Ac 



mi 



16vr 2 M L M E 



Re(\ LE \LE clc R ) ~ c x 10 9 Re(X LE X LE ) 



300 GeV \ 2 mj 
y/M L M E ) mf 



(68) 



where c ~ 0(1). From eq. (68) we learn that if the heavy leptons have masses around the EW 
scale Ml,M e ~ v and if Re(XLE^LE) ~ 1 then the muon g— 2 anomaly can be solved. We 
stress that, in general, the parameters Xle, ^le, Ml^ e , and cl :R are not flavor universal and 
therefore naive scaling might be violated. 

The Yukawa couplings Xle and Xle could also violate lepton flavor and CP. As a reference 
framework, we consider the so-called anarchic scenario, where all the entries in Xle and Xle 
are assumed to be of order one. For the branching ratio of /i — > cy we obtain 



BR(/i -> erf) « c 2 x 10~ 6 | A L£ A LE | : 



/ 300GeV 
WM L M E 




Ml +14/41 



m e /m^ 



(69) 



and therefore, in order to satisfy the experimental constraint BR(^ — > e"f) < 2 x 10 12 , we need 
\/s l t (R A 2 < 2 x lO -6 ^ for c ~ 1. Such a suppression for the flavor mixing angles could be 



achieved by introducing a flavor symmetry determining the structures of Xle and Xle- 
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For the electron EDM we find 

de r 05 t / \ "\ ,/ 300GeV\ 2 

« 5c x 10" 25 Im A L£ A L£ -==== ecm, 70 
e \yjM L M E J 

and the experimental bound is satisfied for \1to.{\le^le)\ ^ 2 x 10~ 3 for c ~ 1. Thus, we 
conclude that both the flavor and CP structures of \le and \le have to be highly non-generic 
to satisfy the constraints from £j — > 1^ and the electron EDM. 

Finally, we mention that interactions with heavy leptons generate also corrections to the 
fermion couplings of the Z and Higgs bosons. Indeed, after integrating out the heavy states 
using their equations of motion we obtain dimension-six operators leading to corrections of the 
fermion couplings of the form v 2 /M L E . 

7 Conclusions 

The aim of this paper was to show that the anomalous magnetic moment of the electron a e can 
be viewed today as a new player among the low-energy processes that are able to probe new- 
physics effects. This novel status of a e stems from recent improvements on both the experimental 
and theoretical fronts. One important ingredient is the measurement of a from atomic-physics 
experiments, which are becoming competitive with a e in the determination of the fine-structure 
constant. The second ingredient is the ongoing effort to measure a e with better experimental 
accuracy. The third element is a more precise theoretical determination of a e in the SM. At 
present, the experimental measurement of a e is in good agreement with the SM and the uncer- 
tainty in the quantity Ao e = «e XP — a^ M is about 8 x 10~ 13 . As discussed in this paper, future 
progress can reduce this error by about one order of magnitude. 

From the theoretical point of view, the great interest in testing new-physics effects in a e 
comes from the well-known discrepancy between the experimental measurement and the SM 
prediction of a„. Observing or excluding an anomaly in a e could become the most convincing 
way to establish the origin of the discrepancy. In a large class of models, new-physics 
contributions to ai (for i = e,fi,r) are proportional to mj, a situation that we call "naive 
scaling". In the case of naive scaling, the present value of the anomaly, see eq. ([I]), 

corresponds to Aa e = (0.7 ± 0.2) x 10~ 13 . While theoretical predictions have recently achieved 
an impressive precision 0(1O -13 ) via formidable calculations, confirmation of the discrepancy 
through a e still requires experimental improvements at the utmost level. 

In many well- motivated cases, new physics effects do not respect naive scaling, as we have 
shown with several examples of such theories. In the context of supersymmetry, this can happen 
with lepton flavor conservation (for non-degenerate but aligned sleptons) or with lepton flavor 
violation (with relative misalignment between lepton and sleptons). In the first case, once we 
normalize Aa e in such a way to reproduce Aa^ ~ 3 x 10 -9 , we obtain specific predictions on 
violation of lepton universality in P — > [iv/P — > ev and P — > fii>/r —> Pv, for P = tt, K. The 
case of lepton flavor violation cannot explain the anomaly in a«, but can lead to new effects in 
a e that are correlated with the prediction of r — > e"f and violation of e/ jj, lepton universality. 
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In a fairly model-independent way, it is also possible to correlate new effects in Aa e with the 
corresponding electric dipole moment. 

We have also considered a class of theories with a light (pseudo) scalar interacting with 
matter proportionally to ordinary Yukawa couplings. In this case we have found an interesting 
pattern of violations of naive scaling coming from the interplay between one-loop contributions 
with A<2£ oc m\ and two- loop contributions with Aap oc mj. For the parameters capable of 
explaining the (g— 2) M anomaly, we found that Aa e is larger than its naive-scaling value and can 
be close to its present experimental bound. 

We believe that a e offers a special opportunity to test new-physics effects and to shed 
new light on the current discrepancy in the muon anomalous magnetic moment. A robust and 
ambitious experimental program is necessary to exploit its full potential. 
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Appendix 

The loop functions entering the supersymmetric contributions Aa^ FC , Aa^ FV , and BR(^j — > £jj) 
are given by 

l-x 2 + 2xlogx 
Ux) = {l _ x y ' ( 71 ) 

_ 3-4x + x 2 + 21ogx 
/cW " (X -1)3 ' [U) 



1 + 4x - 5x 2 + 2x(x + 2) log x 
4(1 -x) 4 



9n (x) = ^ ,,7 T T ^ , (73) 



5-4x-x 2 + 2(2x + l)logx 
Ux) = , (74) 

l + 9x-9x 2 -x 3 + 6x(x + l)logx 
K{x) = , (75) 

, , , 3 + 44x - 36x 2 - 12x 3 + x 4 + 12x(3x + 2) log x 

Ux) = 6(l3x) e , (76) 



We have also defined / (c , n )(x,y) = /( c ,„)(x) - / (c , n )(y). 
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